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ABSTRACT 

We consider the orbital evolution of the S-stars, the young main-sequence stars near the supermassive 
black hole (SBH) at the Galactic center (GC), and put constraints on competing models for their origin. 
Our analysis includes for the first time the joint effects of Newtonian and relativistic perturbations 
to the motion, including the dragging of inertial frames by a spinning SBH as well as torques due 
to finite- N asymmetries in the field-star distribution (resonant relaxation, RR). The evolution of the 
S-star orbits is strongly influenced by the Schwarzschild barrier (SB), the locus in the (E,L) plane 
where RR is ineffective at driving orbits to higher eccentricities. Formation models that invoke tidal 
disruption of binary stars by the SBH tend to place stars below (i.e., at higher eccentricities than) 
the SB; some stars remain below the barrier, but most stars are able to penetrate it, after which they 
are subject to RR and achieve a nearly thermal distribution of eccentricities. This process requires 
roughly 50 Myr in nuclear models with relaxed stellar cusps, or > 10 Myr, regardless of the initial 
distribution of eccentricities, in nuclear models that include a dense cluster of 10 M Q black holes. We 
find a probability of < 1 % for any S-star to be tidally disrupted by the SBH over its lifetime. 
Subject headings: black hole physics-Galaxy:center-Galaxy:kinematics and dynamics-stellar dynamics 



1. INTRODUCTION 

Observations of the Galactic center (GC) reveal a clus- 
ter of about 20 stars, mainly main-sequence B stars, 
that extends outward to about a tenth of a parsec from 
the c entral supermassive black hole (SBH; IGhez et al. I 
2008). These stars, usually referred to as "S-stars," fol- 
low orbits that are randomly oriented and have a nearly 
"thermal" distribution of eccentricities, N(e)de ~ ede 
(jGillessen et al. I I2009D . The existence of such young 
stars so close to the GC SBH challenges our understand- 
ing of star formation since the strong tidal field of the 
SBH should inhib it the collaps e and fragmentation of 
molecular clouds ((Morris] [1993) . For this reason, it is 
usually assumed that the S-stars formed elsewhere and 
migrated to their current locations. However, the mi- 
gration mechanisms proposed in the literature result in 
orbital distributions that differ substantially from what 
is observed. Post-migration dynamical evolution due to 
gravitational interactions with other stars or stellar black 
holes (BHs) has been invoked to bring the predicted or- 
bital distributions more in line with observations (e.g., 
Merritt et aT~l 120091: IPerets et al. Il2009l iMadigan et al. I 



initial conditions that correspond to the different forma- 
tion models proposed in the literature. Evolving the ini- 
tial conditions for a time of the order the lifetime of the 
S-stars, and comparing with the observed distribution of 
orbital elements, allows us to place constraints on both 
the parameters of the nuclear cusp and the S-star origin 
models. 

2. GRAVITATONAL ENCOUNTERS NEAR THE SBH 

Timescales of interest are of order 100 Myr, the main- 
sequence lifetime of a B star, or less. Such times are short 
compared with two-body (non-resonant, NR) r elaxation 
times near the center of the Mi lky Way (e.g., iMerritt I 
120101 : I Antonini fc Merritt I [2012( 1. hence we ignore NR 
relaxation in what follows and assume that orbital ener- 
gies, i.e. semi-major axes a, are unchanged once a star 
has been deposited near Sgr A\ 

Resonant relaxation (RR) (|Rauch fc Tremainel 119961 : 
iHopman fc Alexander 1 12009 1 acts to change orbital ec- 
centricities in a time 



2011b IZhang et al. ll2012ft 



The S-stars approach closely enough to Sgr A* that rel- 
ativistic corrections to their equations of motion can be 
important. In this paper, we apply recent insights about 
how relativity interacts with Newtonian (star-star) per- 
turbations near Schwarzschild and Kerr SBHs. Using 
an approximate Hamiltonian formulation that includes a 
post-Newtonian description of the effects of relativity, we 
explore the evolution of the S-star orbits starting from 
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the "incoherent RR time," where L c is the angular mo- 
mentum of a circular orbit having the same semi-major 
axis as the test star and £ co h is the "coherence time," 
defined as the time for a typical field-star to change its 
orbital orientation; the latter is the shortest of the mass 
precession time (due to the distributed mass), the rela- 
tivistic precession time (due to the 1PN corrections to 
the Newtonian equations of motion), and the time for 
RR itself to reorient orbital planes. For instance, in the 
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Figure 1. Location of the Galactic center S-stars on the (a, e) plane, compared with the Schwarzschild barrier (dashed line, equation [3j , 
and the curve along which frame-dragging torques compete with yN torques from the stars (dash-dotted lines, equation \5\ with two 
diffe rent values of the SBH spin x)- The three panels are for three models of the nucleus, as described in the text. The S-star data are 
from [Gillcsscn ct al. (2009). Horizontal tick marks g ive the expected am plitude of eccentricity changes as an orbit precesses in the fixed 
torquing field due to the field stars (equation 41 from Merritt ct al. 2011). 



case that field-star precession is dominated by relativity, 



y RR ~ — ?JL f Ml\ — — — 

vr 2 a \ m ) N{< a) 

1/2 

1.4 x 10 6 ' 



10 mpc 



I Mr. 



N 

103 



(2) 



yr 



where r g = GM m /c 2 1 P is the orbital period of the test 
star, N(< a) is the number of field-stars with semi-major 
axes less than a, m is the mass of the field stars, and mpc 
is milliparsecs; M, = 4 x 10 6 Mq has been assumed. 

RR ceases to be effective at changing the eccentricities 
of stars whose orbits lie below (at higher eccentricities 
than) the "Schwarzschild barrier" (SB), the locus in the 
(a, e) plane where relativistic precession of the test star 
acts in a time shorter than the time for the field-star 
tor ques to change L. T he SB is defined approximately 



torq ues to cnange l. i 
by ([Merritt et al. II2011I ) 



(1 



1/2 
'SB 



M. 



a m N /7V(< a) 



(3) 



Orbits above (at lower e than) the SB evolve in response 
to RR by undergoing a random walk in e. If such an orbit 
"strikes" the SB, it is "reflected" in a time of order the 
coherence time and random- walks again to lower e, in a 
time ~ Trr, before eventually striking the SB again etc. 
Penetration of the SB from above can occur but only on 
a timescale that is longer than both the RR and NRR 
timescales (jMerritt et al. 112011ft . 

If a star should find itself below the SB, torques from 
the field stars are still able to change the orientation 
of its orbital plane ("2d RR") even though changes in 
eccentricity are suppressed. The timescale for changes in 
orientation is 
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where again M, = 4 x 1O 6 M has been assumed. How- 
ever, 2dRR itself ceases to be effective for orbits that 
come sufficiently close to the SBH, where dragging of in- 
ertial frames by a spinning SBH induces Lense-Thirring 
precession with a period that is shorter than the time for 
2dRR to randomize orbital p lanes. The condition for a n 
orbit to be in this regime is (IMerritt fc Vasiliev ll20T2) 
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with x = cS/(GM 2 ) the dimensionless spin of the SBH, 
and S the SBH spin angular momentum. We define ok, 
the "radius of rotational influence" of the SBH, as the 
value of a that satisfies equation ([5]) with e = 1; ok 
is roughly 1 mpc fo r the Milky Way assuming \ = 1 
(jMerritt et al. Il2010j ). 

While the joint evolution of an ensemble of stars 
near a spinning SBH can only be convincingly treated 
using an iV-body code, Monte-Carlo algorithms have 
been constructed that faithfully reproduce the eccentric- 
ity evolution of single (test) stars due to the dynami- 
cal mechanisms described above, assuming that N(a, e) 
for the field-star distribution is not evolving. In this 
pa per, we use an algori thm similar to that described 
by IMerritt et ai~l (|201lD . The Hamiltonian that de- 
fines the test-star motion includes terms representing 
the effects of the spherically-distributed mass (which re- 
sults in precession of the argument of periastron), 1PN 
(Schwarzschild), 1.5PN (Lense-Thirring) precession due 
to relativity and dipole and quadruple order terms rep- 
resenting the torquing due to the finite-TV asymmetry 
in the field-star distribution (which induces changes in 
all the orbital elements). The direction of the torquing 
field is changed smoothly with time and is randomized 
in a t ime of i co h, as described in §VB of IMerritt et al. I 
(|20ll . 

Above the SB, where test-star precession times are 
comparable with typical field-star precession times, the 
assumptions underlying the derivation of resonant relax- 
ation are satisfied and the algorithm correctly reproduces 
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Table 1 

Origin models for the S-stars 

p" TDs"(%) 



Binary Disruption 



Burst sccnario c 


5Myr 






20Myr 






50Myr 






lOOMyr 




200Myr 




7 = 0.5; m = 1M 
7 = 7/4; m = IMq 
7 = 2; m = 1OM 


7.41 x 
5.11 x 
0.593 


10" 
10" 


-14 
-5 


7.39 x 

0.176 

0.296 


10" 


-12 


2.01 x 

0.765 

0.239 


10" 


IS 


1.31 x 10~ 
7.93 x 10" 
0.160 


s 

2 


9.36 x 10~ u 
3.10 x 10~ 2 
0.201 




1.9 
0.36 


Migration from 
gaseous disk c 


5Myr 






20Myr 






50Myr 






lOOMyr 




200Myr 




7 = 0.5; m = IMq 
7 = 7/4; m = IMq 
7 = 2; m = 10M o 


2.90 x 
3.12 x 
6.78 X 


10" 
10" 
10" 


-11 

-9 
-3 


9.57 X 
6.14 x 
3.76 X 


10" 
10" 
10" 


-10 

-7 

-2 


8.29 x 
4.06 x 
0.138 


10" 
10" 


a 

5 


3.62 x 10" 
1.70 x 10~ 
0.168 


V 

-4 


1.12 X 10~ 5 
1.37 x 10~ 3 
0.184 




0.12 
0.41 




Binary Disruption 
Continuous scenario^ 


5Myr 






20Myr 






50Myr 






lOOMyr 




200Myr 




7 = 0.5; m = 1M 
7 = 7/4; m = IMq 
7 = 2; m = 1OM 


7.21 X 
3.14 x 
8.07 x 


10" 
10" 
10" 


-13 

-5 

-2 


6.41 X 

0.147 

0.108 


10" 


-13 


5.32 x 

0.645 

0.410 


10" 


13 


1.33 x 10~ 

0.819 

0.499 


13 


1.267 X 10~ la 

0.660 

0.310 




0.16 




"Probability value of the 2 samples Kolmogorv-Smirnov test 
6 Percentage of stellar tidal disruptions after 200 Myr 
c Orbits initialized at t=0 

d Orbits initialized at random times between [0,200 Myr] 
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Figure 2. Initial (red points) and final (black open circles, after 200 Myr) locations of of the test particles in the Monte-Carlo integrations. 
Dashed and dot-dashed lines are defined in Figure [T] dotted lines give the tidal disruption radius of a 10 Mq star. Stars that initially have 
large eccentricities and lie near, but to the left of, the SB can penetrate the barrier and move to the right, where they remain. These stars 
end up with a nearly uniform distribution of angular momenta. Stars that are initially right of the SB tend to remain there, though some 
barrier penetration (from right to left) is observed near a cr i t , the limiting value of a for which the SB exists. 
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the eccentricity evolution predicted by RR, as well as the 
"bounce" observed in TV-body integrations when an orbit 
strikes the SB. For a test star that finds itself below the 
SB, the Schwarzschild precession time is short compared 
with field-star precession times. In this regime, a test 
star precesses with period close to 



P a , 

*OR=-r — (1 
o r„ 
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the 1PN apsidal precession time. During one preces- 
sional period, the field-star torques are nearly constant; 
as a result, the test-star's angular momentum oscillates 
with period £gr and with approximate amplitude (in the 
small-£ limit) 



A£fa2eLAu sini, 



A D = 
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M*(< a) a 



3^N{< a) M. 
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(|Merritt et al. 1120111 ). Here, I 2 = 1 - e 2 , M*(< a) is the 
mass in stars at radii r < a, and sini specifies the incli- 
nation of the major axis of the torquing potential with 
respect to the orbit. By themselves, these periodic varia- 
tions in e do not imply any directed evolution in angular 
momentum, but random switching of the direction of the 
torquing potential does result in a random walk in a test 
star's angular momentum, allowing a star that is initially 
below the SB to approach it. 

Dragging of inertial frames results in orbit-averaged 
rates of change of the argument of periastron, w, and the 
angle of nodes, fi, of the test star according to: 



dt J FD 
dt i FD 



2G 2 M 2 
c 3 a 3 (l — e 2 ) 
6G 2 M 2 
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In equations ©, the "reference plane" for f2 and for i 
(the orbital inclination) is the SBH equatorial plane. 

Since little is known about the distribution of stars and 
stellar remnants near the Galactic center, we explored 
a range of different models for the field-star distribu- 
tion. Assuming a power-law density profile, the number 
of stars at radii less than r is 



N{< r) = /Vo.2 



0.2 pc 



3-7 



(9) 



where A .2 = N(< 0.2 pc). The parameters {/Vo.2, 7, m} 
then uniquely define the background distribution in 
which the test particle orbits are evolved. In two mod- 
els, we set m = 1 M©, and we ta ke either 7 = .5 and 
N0.2 = 8 x 10 4 (equation (5.247) inlMerritt 1 12013?) . simi- 
lar to what is inferred from observations (e.g.. IDo et al. I 
|20Q9f). or -v = 7/ 4 and 7V . 2 = 1.6 x 10 5 (equation (5.246) 



in lMerritt1l2013t ) , the expect ed values for a dynamic ally- 
relaxed population of stars (jBahcall fc Wolf I Il977t ). In 
another model, we adopt A .2 = 4.8 x 10 3 , 7 = 2 and 
m = IOMq. This latter choice of parameters approxi- 
mately reproduces the density of stellar BHs predicted 
by collisionally relaxed models of a cusp of stars and 
stella r remnants around SgrA* (Hopman & Alex ander! 
12001 . 



Figure [T] plots the S stars on the (a, e) plane, as well 
as the location of the SB; the latter depends on the pa- 
rameters defining the nuclear cusp through equation ([3]) . 
Dot-dashed lines in the figure delineate the region where 
frame-dragging torques from a spinning SBH would dom- 
inate stellar torques, equation (Q}. Particularly for large 
X, many of the S stars lie close to this transition re- 
gion, suggesting that frame dragging could be an impor- 
tant influence on their orbital evolution; for instance, by 
inhibiting 2dRR. This figure does not give information 
about timescales, but we note that characteristic times 
like Trr are functions of the nuclear parameters, which 
can be important given the limited lifetimes of the S 
stars. 

In the two models with a steep cusp, the SB delineates 
the boundary of the S-star orbits, with only a few stars 
lying below the minimum a = a cr j t ~ (M m /mvN)r g for 
which the barrier exists. In these models, the existence of 
the SB is expected to strongly influence the eccentricity 
evolution of at least some of the S stars; for instance, by 
limiting the maximum eccentricity attainable by a star 
that starts above the barrier. Furthermore, for some of 
the nuclear models, Figure [1] shows that some of the S- 
stars can lie both above a cr ;t and below the SB. Such a 
location would be highly unlikely, in a time as short as 
~ 100 Myr, for stars that started above the SB, but is 
reasonable if the stars were placed initially on such orbits 
via one of the mechanisms described below. 

3. FORMATION MODELS 
We considered two models for the origin of the S stars. 

(1) Formation of the S-stars in binaries far from the 
center (r > 0.1 pc). In this model, the bina- 
ries are scattered onto low-angular-momentum or- 
bits that bring them close enough to the SBH 
that an exchange interaction can occur, leavi ng one 
star on a tightly-bound orbit around SgrA* J Hills 
GUI \Yu fc Tremaine1[2003l: lAntonini et al.ll20ToF 7 
The radius at which the SBH tidally disrupt a bi- 
nary is typically a few tens of AU for main se- 
quence binaries, and the orbital eccentricity is ex- 
pected to be large. The initial orbital inclinations 
of the S-stars will be either randomly distributed 
if the b inaries originated in an isotropic stel- 
lar cu sp (|Perets et al. 1120071 iPerets fc Gualandris I 
2010) or highly correlated if they formed in a stellar 
disk (jMadigan et al. II2009D . 

(2) Formation of the S-stars in a disk at roughly their 
current radius, either one of the known stellar 
disks, or a pre-existing one. Formation in the disk 
would be followed by migration to their current lo- 
cations. This model pred icts initially small eccen- 
tricities and inclinations (|Levin II2007D . 
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Figure 3. Cumulative distribution of eccentricities for different S-star formation models after 5, 20, 50, 100 and 200Myr (line thickness 
increases with time). Open circles give the observed distribution; dot-red line shows a "thermal" eccentricity distribution. Left and right 
panels correspond to formation through capture following binary disruption, middle panels to formation in a stellar disk followed by 
migration through interaction with a gaseous disk. 



4. ORBITAL EVOLUTION 

We start by assuming that N(a) is known: it is given 
by the observed values of a. For each S-star (i.e. for 
each value of a), 100 Monte-Carlo experiments were car- 
ried out using the Hamiltonian model described above, 
in each of the three nuclear models, for an integration 
time of 200 Myr. The initial orbital eccentricities were 
assigned from a thermal distribution, iV(< e) oc e 2 , 
over some specified range in e. In the case of migra- 
tion from a gaseous disk we required e < 0.5 initially. 
When considering the binary disruption model, we only 
considered orbits with initial eccentricities in the range 
0.93 < e < 0.99. In this case, we assume either that 
the S-stars are brought to their current location at the 
same time (burst scenario) or that they arrive at random 
times between and 200 Myr (continuous scenario). The 
former choice corresponds t o a burst of S-star form ation, 
for instance in a stellar disk (|Madigan et al. 1120091 ) . while 
the latter choice assumes that the S-stars form continu - 
ously in the isotropic stellar cusp (jPerets et al. 1 120071 ) . 



In all cases we set \ = 1- 

Stars were assumed to be tidally disrupted when 
they approached the SBH within a distance r t — 

2R (M./m) 1/3 ijAntonini et al. 11201 lal ). with m = 10 M Q 
and R = 8 R Q . 

Figure [5] compares the initial and final (after 200 Myr) 
a — e distributions. In the two, steep-cusp models with 
high initial e, most of the stars start off to the left of the 
SB, where evolution in angular momentum is strongly 
suppressed by the Schwarzschild precession. Neverthe- 
less, it can be seen that, after some time, an initially 
eccentric population separates into two subpopulations: 
stars that are so far leftward of the SB initially that their 
eccentricities hardly evolve; and stars that either begin 
rightward of the SB, and remain there, or that are close 
enough initially to the barrier to cross it. The latter stars 
are subject to RR after crossing the barrier and end up 
with a nearly thermal eccentricity distribution. A clear 
"gap" between the two populations is evident in several 
of the frames of Figure [2j the gap extends from the SB 
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on the right, to a somewhat higher eccentricity on the 
left. The SB acts like a membrane that is permeable in 
one direction only, from left to right: having crossed the 
barrier from left to right, a star moves quickly (in a time 
of ~ Trr ) to a region of lower e where it remains Q. The 
observed evolution below the SB is due to reorientation 
of the torquing potential which results in a random walk 
in angu lar momentum, prefe rentially toward lower e (see 
§VB of iMerritt et al. 1 1201 ll ): this mechanism is qualita- 
tively similar to resonant relaxation but obey s a different 
set of relations (jAlexander fe Merritt"ll2012f) . 

On the other hand, when eccentricities are initially low, 
as in the bottom panels of Figure EH they tend to remai n 
low, i.e., above the SB. As noted iii lMerritt et al. I (| 2 1 If ) . 
the permeability of the SB tends to increase near a C rit 
and this can be seen in Figure E3 as well. 

Time evolution of the cumulative distribution of ec- 
centricities is shown in Figure [3] (In this figure and 
in the analysis that follows, we do not include S- 
stars that likely belong to the disk(s) of O/WR stars, 
i.e. S66, S67, S83, S87, S96, S97). Eccentricity distribu- 
tions were found to approach a nearly "thermal" form in 
a time of order Trr. Based on Figure [3l we see that the 
lifetime of the S-stars may, or may not, be long enough 
for this to happen, depending on the nuclear model. We 
compared the results of the Monte-Carlo integrations 
with S-star data by performing 2-sample Kolmogorov- 
Simornov (K-S) tests on the e distributions (Table [T]). 
The best match to observations is attained after 20 Myr 
of evolution in stellar cusp models in the continuous sce- 
nario with 7 = 7/4 and starting from initially high ec- 
centricities (K-S test p-values of « 0.7). Integrations 
that include stellar BHs also generate orbital distribu- 
tions which are in agreement with observations after ap- 
proximately 5 and 20 Myr of evolution for initially high 
and low eccentricity distributions respectively. 

We tested the degree of randomne ss of the orbita l 
planes using the Rayleigh statistic TZ (|Ravleigh I I191QT ) . 
defined as the resultant of the unit vectors U, i = 
l...iV mc , where U is perpendicular to the orbital plane 
of the ith star and N mc is the total number of Monte- 
Carlo data points (i.e., N mc = 1900). Since the test stars 
were initialized with the same inclination, the orbits are 
initially strongly correlated and TZ ~ N mc ; over a time 
of order T^dRR, 2dRR ten ds to randomize the orbital 
planes and TZ approaches \J iV mc , the value expected for 
a random distribution. When the main contribution to 
dynamical relaxation comes from stellar BHs, TZ reached 
values consistent with isotropy after ~ 20 and ~ 50 Myr 
for initially high and low eccentricity distributions re- 
spectively. In the stellar-cusp model with 7 = 7/4 and 
staring from initial high e, TZ reached values consistent 
with isotropy after ~ 100 Myr, while in all the other 
stellar-cusp models we measured a significant departure 
from randomness (TZ/N mc > 0.2) even after 200 Myr 
of evolution. We estimated TZ separately for orbits that 
any time were below the dot-dashed lines of Figure [2] and 
found that this population had a distribution of orbital 
planes which was less consistent with being random. Evi- 
dently, 2dRR was somewhat inhibited by frame dragging 
for orbits with initially large eccentricities. 

1 This behavior has been confirmed via direct TV-body integra- 
tions. 



The fraction of stars that would have been tidally dis- 
rupted after 200 Myr of evolution was never larger than 
~ 1% (Table [p. As a compa r ison, in purely Newto- 
nian integrations, iPerets et all (|2009f) found that up to 
~ 30 % of stars were disrupted after 20 Myr of evolution. 

5. CONCLUSIONS 

In this paper we studied the combined effects of New- 
tonian and relativistic perturbations on the angular mo- 
mentum evolution of the Galactic center (GC) S-stars. 
For the first time we have shown that the a — e distri- 
bution of the S-stars predicted by the binary disruption 
model, in which the stars are delivered to the GC on 
high-eccentricity orbits, is consistent with the observed 
orbits even when relativistic effects are considered. In 
these formation models, most of the orbits lie initially 
below the Schwarzschild barrier, the locus in the (a,e) 
plane where resonant relaxation is ineffective at chang- 
ing eccentricities. Contrary to this basic prediction, we 
found that orbits starting sufficiently close to the barrier 
are sometimes able to penetrate it, diffusing above and 
reaching a nearly thermal e-distribution; a small fraction 
of stars remain confined below the barrier at low angular 
momenta (e > 0.95). A good match to observations is 
achieved after ~ 20Myr of evolution if the distribution of 
field stars at the GC follows a nearly dynamically relaxed 
form. Models that include a mass-segregated population 
of stellar BHs also generate after ~ 10 Myr of evolution 
distributions that are marginally consistent with obser- 
vations. 

Based on the origin models considered here, the 
S-star obits can only be reproduced by postu- 
lating dynamically relaxed states (i.e., steep den- 
sity cusps) for the GC. This result is interesting 
given that s uch models are curren t ly disfavo r ed by 
observations (jBuchholz et al. I 120091: iDo et al. I 12009 ) 
and by some theoretical arguments (|Merritt I I201C ; 
lAntonini et"aHl2011bl fGualandris fc Merritt 112013 ). 

DM was supported in part by the National Science 
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